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Abstract. A uniform, flexible and frictionless chain falling link by link from 
a heap by the edge of a table falls with an acceleration g/3 if the motion is 
nonconservative, but g/2 if the motion is conservative, g being the acceleration 
due to gravity. Unable to construct such a falling chain, we use instead higher- 
dimensional versions of it. A home camcorder is used to measure the fall of a three- 
dimensional version called an xj/z-slider . After frictional effects are corrected for, 
its vertical falling acceleration is found to be a x /g = 0.328 ± 0.004. This result 
agrees with the theoretical value of a x /g = 1/3 for an ideal energy-conserving 
xj/z-slider. 



PACS numbers: 01.55.+b,01.65.+g,01.50.Zv 
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Figure 1. (a) The falling chain called the xj/-slider that has one end fixed at 
the origin O of the coordinate system. The chain slides along the table edge after 
release. The falling chain is called an xyz-slider when the stationary segment on 
the table parallel to the edge is set back from the edge by a distance z = d. (b) 
The forces acting on the chain bend of the xj/-slider at a sharp table edge. 

1. Introduction 

A uniform, flexible and frictionless chain falling link by link from a heap by the 
edge of a table is a popular example or problem that appears in many textbooks 
of classical mechanics pQ. The solution given in these books is almost always the 
energy-nonconserving one where the chain falls with an acceleration of a = g/3, g 
being the acceleration due to gravity. This solution was given as early as 1857 by 
Cayley [2j El [1] . Recently, Villarino [4] has called our attention to a footnote in 
Routh's treatise on dynamics [5] that refers to problems of this type involving systems 
with variable masses being solved in the private lecture room of the famed Cambridge 
tutor Hopkins as long ago as 1850. The footnote also mentions a solution to our 
falling-chain problem published in 1856 in a treatise on dynamics by Tait and Steele 
[B]. For this reason, we shall call the chain falling from a heap at the table edge the 
Hopkins-Tait-Steele-Cayley (HTSC) falling chain. 

Mikhailov 7 has given a history of dynamics problems involving variable masses, 
especially those that like the HTSC falling chain were made up for or studied by 
undergraduates in Cambridge in the second half of the 19th century in connection 
with the Cambridge Mathematical Tripos Examination [5]. It thus seems useful to 
extend the review given in Wong and Yasui [1] on the HTSC chain in order to clarify 
the historical circumstances under which these dynamics problems were studied in 
Cambridge. This extended review will be given in section [5] 

Wong and Yasui [T] have taken issue with the energy-nonconserving solution of 
g/3 for the acceleration of the HTSC chain. They point out that the full Lagrange 
equation of motion contains an extra term not given by the solution of Cayley and of 
Tait and Steele. This missing term gives the chain tension with which the falling chain 
is pulling on the falling link. When this term is included, the solution becomes energy- 
conserving and the falling acceleration is increased to g/2. The difference in physics 
between the energy-conserving and nonconserving solutions will also be described in 
section [5J 

It is thus of interest to determine experimentally which of these two solutions for 
the HTSC falling chain is correct. We immediately run into the difficulty that a heap 
or coil of chain on a table cannot be built in a reproducible way, for every heap seems 
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unique. After much thought, we discover a two-dimensional realization of the falling 
chain called the xy-slider that is easily reproducible. This is a chain stretched out 
on the table parallel to the table edge and very close to it, as shown schematically in 
figure []J a). One end of the chain is held fixed on the table with the fingers of one hand 
while the other end overhangs the table slightly. When the overhanging end is allowed 
to fall, the chain is found to slide along the slightly round table edge to the fixed end, 
as shown in the figure. The reason for holding down one end is to ensure that the 
resting 'heap' on the table (now represented by a straight chain segment parallel to 
the table edge) remains stationary on the table while the rest of the chain falls. It also 
makes the chain less susceptible to perturbations caused by the release of the falling 
end. The chain motion then becomes more easily reproducible. 

The motion of a uniform, flexible and frictionless xy-slider is described in section[3j 
using Lagrangian mechanics. Such an ideal chain experiences no dissipation so that 
its energy is conserved during the chain fall. For the special case where the fallen 
chain segment is a straight line at a constant angle 9 to the horizontal table edge, the 
theoretical falling acceleration a is found to be constant in time, with 



where s is the length of the fallen chain segment. 

We now need to determine experimentally how an actual xy-slider falls. In the 
rest of the paper, we shall show that the required measurements can be made using 
only simple equipment found at home, so that the interested reader can readily repeat 
our measurement. We shall begin, in section [5J by showing that an inexpensive home 
camcorder [9] recording at 30 frames a second has sufficient resolution to time the fall 
of chains. We verify its capability by measuring the known value of g for the free fall 
of a metal object. The measurement on the xy-sli&er is finally described in section [5] 
We explain why a three-dimensional version of the slider called an xyz-slider is used. 
The results roughly corrected for friction in the air, on the table, and at the table edge 
are 



The corrected motion of a real xyz-sli&ei is thus found to agree with that of the ideal 
energy-conserving rry-slider. Our result thus gives indirect experimental support that 
the ideal HTSC chain is also energy conserving, and can be expected to have a falling 
acceleration of g/2. 

2. The history and physics of the HTSC falling chain 

According to Mikhailov 7 , rocket propulsion by the continuous emission of masses was 
studied as early as 1810 by Moore. A Lagrange equation of motion for a more general 
mechanical system with variable mass was given by Buquoy in 1812, and re-written in 
a more modern form by Poisson in 1819. We are interested here in the variable-mass 
problem of a flexible and frictionless chain falling link by link from a heap at the table 
edge. This problem is a special case of the following problem described and solved by 



9 = cos" 1 ^) = 70.5°, 
a =s = g/V8 = 0.354 ff , 
a x = asin# = g/3, 



(1) 



9 = 70° ± 2°, 

a x /g = a$in9/g = 0.328 ±0.004, 

a/g = 0.350 ±0.006. 



(2) 
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Figure 2. The falling chain of Tait and Steele \§\. 



Tait and Steele in a textbook on classical dynamics published in 1856 [6, 8, 4 : 'One 
end, B, of a uniform heavy chain hangs over a small smooth pulley A, and the other 
end is coiled up on a table at C. If B preponderates, determine the motion.' This 
falling chain is shown in figure [2j 

Tait and Steele solve the problem by using Newton's force law written in the form 



— [p(x + h)v] =pg(x-h), 



(3) 



where p is the uniform linear mass density of the chain, x is the length AB, the chain 
end as measured from the pulley, h is the length AC, the height of the pulley above 
the table, and v = dx/dt. By multiplying both sides of the equation by (x + h)v, ^ 
can be integrated immediately to give 

\{x + h) 2 v 2 = § - x )(x 2 + x Q x + x 2 - 3ft, 2 ), (4) 

if initially x — xq and v = 0. If xq = 2h also holds, the acceleration becomes constant 
in time, a — g/3. Our falling HTSC chain corresponds to the special case h = xq = 
of this problem, the case of a bare chain without the pulley. 

If the total mechanical energy of the falling chain is conserved, however, the result 
would have been v 2 = gx, with a — g/2. Sommerfeld in his textbook on mechanics |10j 
calls the missing mechanical energy AE = px 2 g/6 Carnot's energy loss to emphasize 
the picture that the missing energy is lost when a falling link from the table falls by 
sticking onto the falling chain in an inelastic collision. Inelastic impacts where the 
masses stick together after the impact are known to lose total mechanical energy, a 
result first described by Lazare Carnot in a work published in 1783 [TUl ITTl fT2] . 

A year after the publication of the treatise of Tait and Steele, Cayley wrote a 
paper [5] in which the same equation of motion is obtained from the Lagrangian 



L\x, v) — -xv H x 

v ; 2 2 



(5) 



We are unable to follow Cayley's arguments, and shall just describe his final result. 
His version of Lagrange's equation contains the same term describing the rate of 
momentum change as the left-hand side of (O, but the part where the variable mass 
of the falling chain is differentiated is obtained in a different way with the help of 
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another function (denoted K in his paper). His prescription is such that in calculating 
the generalized force 

dC(x,v) 

— g- — = pgx + T, (6) 

the x dependence of the variable mass of the first (or kinetic-energy) term of the 
Lagrangian is not to be differentiated. The consequence is that the Wong-Yasui 
tension T = pv 2 /2 that comes from such a differentiation does not appear. As a result, 
Cayley's equation agrees with © of Tait and Steele. Cayley nevertheless considers his 
method more general. Mikhailov [7] has pointed out, however, that Cayley's unusual 
prescription of selective differentiation has not been adopted by any subsequent worker 
in the field. 

It is interesting that Cayley did not mention the Tait-Steele solution in his 
paper. Cayley was very knowledgeable, of course, about theoretical developments 
in dynamics. In fact, earlier in 1857 he wrote a rather detailed and very well-received 
42-page report on theoretical dynamics that summarized the progress of the subject 
since Lagrange's Mecanique Analytique of 1788. In this report, Cayley discussed the 
work of Lagrange, Poisson, Hamilton, Jacobi, Liouville, Bertrand, Denkin and others, 
up to early 1857 [T3l [14], all by name. We believe that the probable reason Cayley 
did not refer to Tait and Steele is that both the problem and its solution were so well 
known in Cambridge circles that a reference to Tait and Steele seemed unnecessary, 
and perhaps even unjustified. 

In this connection, it is worth pointing out that according to Routh 5j , problems 
of impact or infinitesimal impulse like that attributed to our falling chain had 
been studied as early as 1850 in the private lecture room of William Hopkins, a 
popular private tutor ('coach' in Cambridge undergraduate slang [15j ) to Cambridge 
undergraduates preparing for the Mathematical Tripos examination. Cayley, Tait and 
Routh were all coached by Hopkins, who helped them win the highest score in the 
Mathematical Tripos of 1842, 1852 and 1854, respectively [16]. (The highest scorer 
in the Mathematical Tripos each year was known for life as the Senior Wrangler, 
Wranglers being holders of first-class Cambridge B.A. degrees in mathematics.) Since 
Routh was a student in Hopkins's lecture room before the Mathematical Tripos Exam 
of 1854, the date 1850 given in his footnote on Hopkins's work on variable-mass 
problems can be considered reliable. 

On the other hand, Cayley took the Tripos exam eight years earlier in 1842. So it 
appears likely that he did not, like Routh and Tait, learn about this class of problems 
with variable masses directly from Hopkins in the latter's lecture room. We should 
add that after a four-year fellowship at the Trinity College, Cayley left Cambridge 
in 1846 to study law in London [T7] [IS]. However, Cayley in his London days was 
in touch with Cambridge. He served as a Cambridge college examiner in the years 
1848-51 [17]. In 1852, he was elected a fellow of the Royal Society with the help of 
many Cambridge supporters, including William Thompson (Lord Kelvin), Stokes and 
Hopkins [18] . He also served as the Senior Examiner of the Mathematical Tripos that 
year. Tait, the Senior Wrangler that year, became one of his friends [T7]. Cayley's 1857 
paper, received by the journal on 18 June, 1857, begins with the sentence 'There are a 
class of dynamical problems which, so far as I am aware, have not been considered in 
a general manner.' His paper was therefore written about a year after the publication 
of the Tait and Steele book on dynamics 6J, where Tait's Preface was dated January, 
1856, and written from Queen's College in Belfast where he was the professor of 
mathematics since 1854. We give all these historical details in order to suggest that 
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Cayley did not acknowledge Tait and Steele perhaps because the problem and the 
solution were not invented by them. We are not sure if this is true, because we have 
not been able to find any earlier reference on this falling-chain problem. 

Wong and Yasui pQ have given a list of other textbooks on mechanics where the 
energy-nonconserving solution for the HTSC falling chain can be found. Since the 
problem is so closely connected to the Cambridge Mathematical Tripos Examination, 
it is interesting to add to this list of books the 1878 second edition of Wolstenholme's 
popular problem collection [19] for the Mathematical Tripos exams, Jeans's textbook 
on theoretical mechanics [5D] and Lamb's book on dynamics [ST]. Wolstcnholmc (Third 
Wrangler of 1850), Lamb (Second Wrangler of 1872 and a student of Routh) and Jeans 
(Third Wrangler of 1898) all took the Mathematical Tripos Examination [22] , 



3. Motion of the xy-slider 

We have noticed that the xy-slider shown in figureflja) falls in a very special way. Soon 
after release, the falling chain assumes an equilibrium shape in which a rather straight 
falling chain segment slides along the table edge at a constant angle 9. Depending on 
the chain and on the manner of release, there are in general some transverse vibrations 
and distortions from the linear shape. We shall limit ourselves to the simple case 
where these complications are either absent or small, and the fallen chain segment 
is approximately a straight line of length s at an angle 9 to a sharp table edge, as 
illustrated in the figure. For a chain of length L, the coordinates of the falling chain 
end as seen from the fixed end O are 

{x,y,z) = (ssin9,L — s + s cos 9,0). (7) 

Using the two-dimensional form (ps/2)(x 2 + y 2 ) of the kinetic energy of the fallen 
chain segment, we find the two-dimensional Lagrangian 

£ 2 (s,v = s) = psv 2 {\ - cos(9) + ys 2 sin(9, (8) 

where p is the uniform linear mass density of the chain, and v = s. There is no 
contribution from the stationary chain segment on the table as its kinetic energy 
vanishes, and its potential energy can be taken to be zero. Our model is thus greatly 
simplified from the general case of two-dimensional motion, for which each term of 
the Lagrangian is an integral in ds. 

The two-dimensional Lagrangian (|8|) is, like its one-dimensional counterpart, not 
an explicit function of the time t. Hence the total mechanical energy of the system is 
conserved. If initially s = and v = 0, then 

2 g sin9 

v = -s. (9) 

2 1 - cos 9 w 

Such a velocity comes from the constant acceleration 

q sin 9 

a=s = - (10) 

4 1-cos6> v ; 

along the longitudinal direction of the falling chain segment. These solutions for v 2 
and a satisfy the Lagrange equation of motion 

(2ss + s 2 ) (1 - cos 9) ~ sgsin9 = 0. (11) 

To understand the preferred direction 9 of the xy-slider, we consider first a plumb 
line whose support moves horizontally with the same acceleration g' = d 2 (L—s)/dt 2 = 
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~a as the bend of the chain at the table edge. In the laboratory frame, this plumb 
line does not point vertically downward, but at an angle 9 to the table edge such that 

a sin 9 

cot 9 



g 4(1 -cos ey 

(3cos0- l)(cos0- 1) = 0. (12) 

The solution cos 9 — 1 or 9 — is unphysical because it requires a — oo . The physical 
solution is 9 = cos" 1 (1/3) = 70.5°. 

It can be expected intuitively that the chain actually falls along the modified 
plumb direction g' found here. To show that this is indeed the case, we concentrate 
on the variable-mass chain segment of length s that has already fallen below the table 
edge. If the chain bend at the table edge is sharp, the fallen chain segment is acted 
on simultaneously by its own downward weight, an upper chain tension T u along the 
negative y-axis, and a lower tension Tg (if any) pointing up at an angle 9 above the 
horizon, as shown in figure []Tb). Its Lagrange equations in the x- and y-directions are 
then 

psg-T e sin 9 = ^-(psi), (13) 
at 

-T u -T e cos 8 = j t (psy). (14) 

They can be simplified to the Newton equations 

psg — Tn sin 9 = ps'x — psa sin 8, (15) 
A — Tn cos 8 = psy = — psa(l — cos9), (16) 

where A = —T u + ps 2 and the lower Newtonian tension 

T N =T e + ps 2 (17) 

is an effective force that differs in general from the Lagrangian tension Tg. These 
Newton equations can be solved for a and TV: 

A 

a — g cot 9 , 

ps 

T N =psg^—^ + A. (18) 
sin 6* 

The magnitude T u of the upper chain tension at the table edge satisfies the Newton 
equation T u — pL\Y\, where Y is the y-coordinate of the center of mass: 

Y = ±-lL 2 -s 2 (l-cos9)}. (19) 

Hence 

T u = 3/ras(l — cos0), 

A = pas(3cos8- 1). (20) 

where use has been made of (|9|). These expressions when inserted into (|18[) give 

9 9 
a = , or a x = — , 

3 sin 9 3 
T N = 2pas, or T t = 0. (21) 

We have thus shown that the Lagrangian chain tension vanishes below the chain 
bend. The equation of motion (|13[) then becomes identical to that for the incorrect 
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energy-nonconserving solution of the one-dimensional HTSC chain. This is why the 
energy-conserving xy-slider has the same downward acceleration g/3 as the incorrect 
solution of the one-dimensional HTSC chain. We recall that the correct result of g/2 
for the energy-conserving HTSC chain is obtained only with the help of a nonzero 
downward-pointing Lagrangian tension that has been called the Wong-Yasui tension 
in section [2] 

The constant acceleration a of the xy-slider must also satisfy (fTUj) because of 
energy conservation: 

q gsin6> . . 

(22) 



3 sin 6 4(1 -cos 0) 

This requirement can be used to show that the angle 9 satisfies the same condition 
(|12| as the plumb line. The chain thus falls along the plumb direction. When this 
happens, the Newtonian chain tension below the sharp bend has the same magnitude 
as the tension T u above the bend, namely T/v = T u , while A vanishes. 



4. Testing the experimental setup for free fall 

We now have to measure experimentally the acceleration of a real xy-slider. The 
falling time over a distance x under a constant acceleration a is 

t=v /|. ,23, 

For x = 0.5 m, the falling time is t = 0.32 (0.45, 0.55) s when a = g (g/2, g/3). For 
a camcorder recording at 30 frames a second, the event lasts 9.5 (13.5, 16.5) frames. 
Such a camcorder thus appears adequate for our timing needs. 

To test the camcorder as a timer, we first measure the falling time of several 
identical 5/16" steel nuts (as in nuts and bolts) over a height of 1.1-1.3 m. A small 
golden plastic bead is tied to each nut to improve its visibility in the video record. 
The nut is initially held between two fingers, one from each hand. At t = Q, the fingers 
are pulled away horizontally from each other. By extrapolating the recorded finger 
separations in two succeeding frames, we can estimate the starting time to about 0.2 
frame, or 0.007 s. A convenient final frame is then chosen and the falling distance 
shown in it is measured. 

The distance measurement can be done by transferring the 2.5" LCD screen image 
to a PC using the camcorder manufacturer's software. The metal nuts are dropped 
just in front of a lunch counter. The known height of the counter from table top to 
the top of the footrest both showing in the picture allows the measurements to be 
converted into actual distances in meters. 

The falling positions are measured from the video record by two persons 
independently. They both report a reading uncertainty of 0.004 m. (After the 
measurements are completed, we find that the PC images can be captured by another 
software and enlarged for more accurate readings.) All results and errors reported 
here are the averages from these two independent readings. 

The observed acceleration calculated from the total elapsed time turns out to be 
larger than expected. To find out the reason, we measure the falling position x(t) 
frame by frame in the video recording. The data are fitted well by the functional form 

x(t) = c + cit + c 2 t 2 . (24) 
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Table 1. The vertical acceleration measured for free fall and for the xyz-slider. 
The position of the fixed end of the xj/2-slider is also noted. The measurements 
on the xyz-slidei are described in section [5] 



Quantity a 


Simple mean a(8a) 


Weighted mean a(5ga, S^a) 


Free fall g 


9.820 (0.034) 


9.818 (0.017, 0.033) 


xj/z-slider a x : 






Left 


3.149 (0.013) 


3.147 (0.013, 0.011) 


Right 


3.136 (0.016) 


3.129 (0.014, 0.014) 


Average 


3.144 (0.010) 


3.139 (0.009) 



The fit gives a small initial downward velocity c\ > in all runs, thereby showing 
that the person releasing the steel nut does not do it cleanly. The acceleration due 
to gravity obtained from g = 2c2 represents the value after correcting for the initial 
velocity. The corrected values are found to be more consistent between runs by a 
factor of 2 than the values obtained from the total elapsed time. 

A simple (unweighted) average over N = 14 runs gives the experimental result of 
g±6g — 9.82 ± 0.03 ms~ 2 shown in table [TJ Here 5g is the standard deviation for the 
mean of N measurements: 

s 2 1 N 

(^f = ^= W ^ T) Y i (9 i - 9 f, (25) 

where s g is the standard deviation for a single measurement. This Sg contains only 
statistical contributions. No systematic errors will be estimated for free fall. The 
camcorder timer is assumed to have the accuracy 6 x 10~ 6 of a quartz clock [23j . 

A more sophisticated estimate of g is possible. The least-square fit to the frame by 
frame data for each free fall gives a fitting uncertainty Sgoi , taken to be the standard 
deviation or 'standard error' calculated by the Mathematica linear-regression function 
Regress [24] . The weighted mean and its raw standard error are then calculated from 
the formula |25j 

g ±S g = ^fi±(J2^)- 1/2 , (26) 

with the weights w% = (Sgoi)~ 2 . The goodness-of-fit parameter for all N 
measurements, 

X 2 = J2 w ^9i-9) 2 , (27) 

i 

turns out to be larger than N — 1, the number of degrees of freedom. This means that 
the fit is poor, or that the estimated errors are too small. Following standard practice 
[25] . all errors are next increased by a scale factor 

1/2 



k = 



X 2 



N - 1 



(28) 



so that the re-computed value of x 2 /{N — 1) becomes exactly 1 after the scaling. In 
our final report, we use the larger of the two errors (5o<? and 

S k g = k5 g (29) 



for any k. 
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Table [T] shows that the weighted mean agrees very well with the simple mean 
in both value and uncertainty. Their close agreement occurs because the values of 
both gi and Sg^i are clustered together, thus indicating that the experimental data 
are consistent from run to run. The measured value of g = 9.82 ± 0.03 ms -2 agrees 
with the known sea-level value at Los Angeles (33°56' N) of g = 9.796 ms -2 [26]. The 
use of our camcorder as a timer is thus validated. 



5. Result for the xyz-slider 

A real falling xy-slider differs in many ways from the theoretical idealization described 
in section [3] First of all, we need to use a rounded table edge to make sure that the 
chain links round it smoothly. The chain must then fall over the rounded table edge 
link by link rather than collectively. Such a controlled fall is achieved by setting the 
chain back from the rounded table edge by a parallel distance d. The setback xy-slider 
is thus an xyz-slider, with the z-axis on the table and perpendicular to the table edge, 
as shown in figure 1(a). Initially the chain runs along the y-axis from the fixed end at 
the origin O. It then turns 90° to run along the z-axis for a distance d before it runs 
over the table edge to overhang the table top a distance Xo along the vertical cc-axis. 
On release, the moving chain rapidly assumes a shape shown in the figure where the 
setback segment makes an angle tp with the table edge, and has length 

s d = ^-. (30) 

smip 

The falling chain segment in the xy-plane has length s and soon makes an angle 9 
with the table edge. The length of the stationary chain segment on the table along 
the y-axis is L — s — Sd, assuming a sharp table edge for simplicity. 

We now show that a setback d of the chain from the table edge can be compensated 
by a proper choice of the initial overhang xq. With a sharp table edge, the chain 
Lagrangian for the xyz-slider is 

£3(3, v) « p(s + Sd)v 2 (l — cos 9) + tj-s 2 sin 9, (31) 

where we have ignored the difference between the angles ip and 9 and have treated the 
segment Sd as a continuation of s in the kinetic energy. Energy conservation can be 
used to find the longitudinal velocity v — s of the fallen chain segment: 

v 2 M 2a S —^, (32) 
s + s d 

where a = (y/4) sin 9/(1 — cos 9) is the constant acceleration without setback or 
overhang. We see that for the choice so = xo / sin# = Sd, the chain falls with the same 
constant acceleration a as the xy-slider without setback. (|32|) can then be written 
more simply as 

x — xo as I (a sin 9)t 2 7 (33) 

assuming that the overhang angle 9 is constant in time. 

The sharp table edge used in our theoretical analysis greatly simplifies all 
formulas, but it is mathematically different from the rounded table edge used in the 
experiment. In the motion of another chain sliding as a whole off the table edge, we 
have studied equations of motion for both sharp and rounded edges and shown that 
while they arc mathematically distinct, they give states of motion that are very close 
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to each other numerically [57] . Hence the use of a sharp table edge in our theoretical 
analysis is acceptable. 

A real xyz-slider also experiences frictional drags as it moves. There are at least 
four sources of friction: (1) between the table surface and the setback chain segment 
of length Sd on the table, (2) between the table edge and the chain rounding it, (3) 
between neighboring chain links at the bends both on the table and at the table edge 
(bending friction) , and (4) air friction resisting the horizontal and vertical motions of 
the falling chain segment. The influence of both setback and frictional forces can in 
principle be compensated by a proper choice of the initial overhang xq of the chain. 
Frictional force # 1 is just the usual 

Fi = fikpgsd, (34) 

where /!& is the coefficient of kinetic friction. The other frictional effects are much 
more complicated. They will be described in | Appendix A[ 

We have used a variety of light ball and cable chains in our measurements, each 
chain about a meter long. We have used a steel desk, a marble counter top, and 
construction lumber with edges routed to different radii of curvature. Our most 
consistent results are those where the falling chain shows the least amount of transverse 
distortion and vibrations. They are obtained with a light steel cable chain used for 
making necklaces and jewellery. A small golden plastic bead (of diameter 3/16" « 0.48 
cm) is tied to the falling end to make it more visible in the camcorder picture. Chain 
and bead together is 0.949 m long, 11.744 g in mass, to which the bead contributes 
0.0569 g. The chain has 316 links, making the average link 0.299 cm long, and 0.0370 
g in mass. The bead is about 1.5 chain links both in mass and in length. 

The table edge used in the reported result is made from a piece of 2" x 6" 
construction lumber (hereafter referred to as 'the table') with the edge routed by 
a 1/4" pa 0.64 cm (radius) router bit. Both the edge and the wood surface on which 
the chain slides are covered by Teflon tapes [28 to reduce friction. The use of wet 
lubricants on our light chains seems to increase friction through surface tension rather 
than to reduce it. The use of dry graphite powder reduces friction marginally, but 
we decide that it is not worth the mess. So all measurements are done on a dry 
unlubricated surface. Our experience is that construction lumber made of fir is too 
soft and grainy to make a good table edge. We recommend using pine or oak instead. 
For a semi-quantitative demonstration, however, the choice of a table edge is not 
critical. 

The levelness of the tabletop is adjusted by shimming with folded pieces of paper 
and checked by using an inexpensive level embedded in a heavy 4-foot plastic yardstick. 
The frictional coefficient fik for our chain and surface is measured by placing the chain 
on the table in a straight line perpendicular to the edge and partially overhanging 
it. The overhanging length is then varied until the whole chain begins to move. The 
resulting static frictional coefficient is /i s = 0.16 without tapping, and 0.15 with gentle 
tapping of the table. The kinetic frictional coefficient can be a complicated function 
of the velocity [29j [30] , but is usually considered to be significantly smaller than /i s 
[50] . We shall use a relatively large value of ^ = 0.13. 

In our measurements, we first arrange the offset segment along the z-direction 
and the overhang segment in the vertical direction. The chain is held at the bend 
on the table with a pen controlled by hand. The chain is released by rapidly moving 
the pen upwards in a vertical direction. The starting time is determined by linear 
extrapolation of the vertical position of the releasing pen in two succeeding video 
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Figure 3. Time dependence of the overhang angle 8(t) in one of the runs of the 
xj/z-slider. The initial value of 9 is 90°. 



frames. 

The theoretical 'edge' of the rounded table edge is defined to be the line parallel 
to the y-axis at an angle of cf>o ~ 23° from the upward vertical (negative cc-axis). 
This angle is obtained as follows: A point mass initially at rest at 4>o moving on a 
frictionless circle that is the rounded edge in cross section will fall off the circle and 
fall down by gravity to the horizontal position of <fi — 7r/2. The angle c/)q is so defined 
that the falling time of the mass is the same as the time of free fall from the table top 
through a vertical distance that is the radius of the circle to the same final height. 

In doing the measurements, we feel most comfortable working with an offset of 
4 chain links, or d ~ 1.2 cm. The initial overhang is thus xq = d = 4 links for the 
initial values of 6 = ip = 90°. The equilibrium shape after release has ip « 45° as 
estimated by eye, and an overhang angle of 8 = 70° ± 2° as measured from the video 
recording. The compensatory overhang at the equilibrium shape would have been 
xq = d sin 9/ sanp sa 6 links. Since we cannot change the overhang once the chain is 
released, we use a middling value of 5 links and add another 0.5 link for friction. After 
the measurements, a more detailed study of frictional effects shows that this frictional 
correction is too small. So an additional correction will be applied to the measured 
result. 

On release, the change in chain shape begins at the bend on the table with the 
angle tp rapidly taking the value of « 45°. This chain direction extends to below the 
chain bend at the table edge before enough of the chain has fallen down to enforce the 
equilibrium value of sa 70°. The time dependence of 8(t) for one of the runs measured 
by one reader is shown in figure [3] We can see an initial transient in overhang angle 
that lasts 0.2 second. The slight modulation of 6(t) after 0.2 s indicates vibration 
caused by edge imperfections or disturbance on release. 

It is hard to ensure that our chain release is clean. Furthermore, the larger static 
friction acting initially may also affect the result. To reduce these uncertainties, we fit 
the falling distance x(t) frame by frame to the three-term power series (|24j) used for 
free fall. This fit corrects specifically for differences in initial position and velocity in 
different runs, while random errors in the fitted acceleration itself will be reduced by 
averaging over several runs. We do 10 runs with the left chain end held fixed on the 
table and 10 runs with the right chain end held fixed, to average out any unlevelness 
of the table. 

A final complication is that the vertical scale of the video image is not exactly 
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the same from left to right. This change of scale is handled by measuring the vertical 
falling position x in each frame relative to the fiducial distance (between the table top 
and the top of the footrest) at the appropriate y-position of the moving chain end seen 
in that frame. The results for a x from the two independent readers agree well, in fact 
better than the results for free fall. Again, all reported values of a x are the averages 
of the results from the two readers. The accelerations obtained by fitting x(t) to (f24f 
are summarized in table [TJ The average of the two sets of runs with different fixed 
end gives a x = 3.139± 0.009 ms" 2 or a x /g = 0.3204± 0.0010, where the uncertainties 
are statistical only. 

We now turn to systematic corrections and errors. The analysis is rather lengthy, 
and shall be given in |Appcndix A| The conclusion is that the main source of energy loss 
is air friction on the sliding chain and that the raw acceleration should be increased 
by 2.5% ± 1.2% to the value of a x = 3.22 ± 0.04 ms~ 2 or a x /g = 0.328 ± 0.004, 
and a/g = 0.350 ± 0.006. We have not included a correction for vibrations and other 
transverse motions of the chain that is likely to bring the result closer to the theoretical 
value of a x = g /?> for the ideal frictionless chain. 

Our error analysis shows that some of the refinements used are not necessary in 
a semi-quantitative demonstration since we have obtained quite good results with the 
many tables used. Further improvement of the quantitative result obtained will require 
a better correction for air friction, a cleaner method of chain release to reduce chain 
distortion and vibrations, a correction for the energy lost to these transverse motions 
and a better table edge. The effect of the initial transient in chain shape should be 
understood. It would be great to automate the reading of the falling position x{t) 
for it is a rather tedious job. However, our error analysis shows that we have enough 
position resolution. 



Appendix A. Systematic errors and corrections 

The length of the moving chain segment on the table is Sd = d = 4 links initially, but 
Sd = d/ simp ss 6 links on reaching equilibrium shape. A middling value of 5 links is 
used for the purpose of data analysis. We estimate that this choice has an uncertainty 
of 0.5 link. The resulting uncertainty on the experimental acceleration a oxp can be 
found by writing ([32]) in the form 

^raw — ^cxp , i Ixx.i. ) 

s + s d 

where a raw is the raw acceleration calculated under the assumption of exact 
compensation between overhang and setback. This equation shows that if the overhang 
so is too small by an amount Asq = Sd — so, then the measured a raw is too small. It 
should be increased by the fractional amount 

a r aw S + S 

Taking s to be half of the average falling distance, or (0.37/ sin 6) m = 0.39 m, we find 
an uncertainty of Aa/a raw = ±0.4%. 

Of the frictional forces acting on the xyz-slider, the simplest is the kinetic friction 
F\ = ^kPgsd on the moving chain segment on the table. Using the equilibrium value 
of Sd = 6 links, we find a result of fXk^d ~ 0.78 link in units of pg. The frictional 
force at the table edge is F2 = fikN, where the normal force N lies in the xz-plane 
for any angle ip. (The fact that N y vanishes for an ideal table edge that is uniform 
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in the y-direction comes about because there is no preference for either the positive 
or the negative y direction.) This normal force is well approximated by a single force 
acting at a sharp table edge at an angle 4>n = 45° from the upward vertical direction. 
In the absence of any friction on the table or at the table edge, N is related to the 
z-component T z of the chain tension just above the table edge as 

Ncoscj) N w T z = pLZ = pd-fj=, (A.3) 

V8 

where Z is the z-coordinate of the center of mass of the chain. Hence N sa pdg/2. 
If table and edge frictions are included in the force balance in the z-direction, the 
resulting expression is considerably more complicated, as we shall show in |Appcndix B| 
The normal force is then increased to Q.85pdg, and yields the table edge frictional 
force of F2 cx 0.85/ifed = 0.44 link worth of chain weight. In the experiment, we have 
allowed only 0.5 link for these frictional forces. There is thus a shortfall of 0.72 link. 
Consequently the raw acceleration a raw should be increased by 0.5%. 

To understand air friction F4. on the sliding chain moving in air, we study its 
effect on the motion of a certain falling chain studied experimentally by Tomaszewski, 
Pieranski and Geminard [3T] (TPG) and shown in their figure 6(a). This chain is a 
light ball chain that is initially stretched taut in the form of a horizontal flat catenary. 
One end is then released while the other end is held fixed. The falling end is found to 
be part of a horizontal chain segment that is falling freely. The length of this freely 
falling horizontal chain segment decreases until it merges completely into the swinging 
arm attached to the fixed support. The almost straight swinging arm at that moment 
has swung down from the horizontal direction by an angle of roughly 80°. This is just 
after frame # 22 in TPG figure 6(a) where the (vertical, horizontal) displacement of 
the chain end from its initial position is (h,w) = (0.96,0.76) m. 

To isolate the effect of air friction on this TPG falling flat chain, we note that 
the time of free fall to frame # 22 is t FF = ^2h/g = 0.442 s. We find from TPG 
figure 5(a) that the their theoretical falling time for that frame is smaller than their 
experimental time by At = 0.0026 ± 0.0003 s. The TPG theory already contains 
bending friction between links in the chain, but not air friction. Hence the fractional 
effect of air friction is 

-^ = 0.0059 = -^, (A.4) 

where a (= g here) is the constant vertical acceleration. (We do not use the free-fall 
curve in TPG figure 5(a) to estimate air friction for the reason that this TPG chain 
actually falls slightly faster than g because its initial position is not exactly horizontal.) 

The TPG falling flat chain is simply related to our xy-slider. The xy-slider also 
has a horizontal chain segment, the part still on the table, that diminishes in length 
as the chain falls. However, the ly-slider is fixed at the end of the horizontal segment, 
while the falling flat chain is fixed at the end of the swinging arm. So they look very 
much like an upside-down version of each other, but with some important differences. 
The angle 9 of the swinging arm in the falling flat chain changes continuously, instead 
of being constant for a significant part of the fall. What is much more important to 
us is that these chains have different falling accelerations: the TPG chain falls with 
a vertical acceleration g, but the xy-slider falls with a vertical acceleration g/3. How 
does this change in acceleration affect the fractional change in the falling time when 
the retarding effect of air friction is removed from the experimental result? 

This question is easily answered for a mass m falling down vertically with 
constant acceleration a = x that experiences an additional weak frictional force 2j3mv 
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proportional to the instantaneous velocity v — x like that acting on the Millikan 
falling oil drop. The resulting acceleration change, Aa — —2f3v, can be expressed as 
the fractional change 

/(a) = 1^1 = 2pt = 2pM. (A.5) 
a V a 

Hence the fractional acceleration changes of the same mass falling through the same 
vertical distance x in the same medium under two different accelerations a\ and a-i 
are related as 

\faif{a-i) = const = y/a^f(ax). (A. 6) 

If the TPG chain falls vertically like this mass, we shall have 

f(g/3) = V3f(g) = 0.020. (A.7) 

This results holds for the xyz-slider also if its factional coefficient j3 is the same as 
that for the TPG chain. 

The TPG chain is a light ball chain, while ours is a cable chain that is shorter 
in length and smaller in mass, 0.95 m and 12 g instead of 1.02 m and 21 g. Their 
cross sectional dimensions are comparable, but their cross sectional shapes are not the 
same. However, both our xyz-sWdei and the TPG flat chain are falling at roughly the 
same angle to the vertical direction when air friction is strongest, i.e., when the falling 
speed is greatest towards the end of the measured fall. So their effective cross sections 
are quite similar in spite of the differences in shapes and in falling configurations. 
We conclude that their (3 coefficients are likely to be quite similar, and that our 
experimental vertical acceleration should be increased by 2.0% before comparison 
with the theoretical value for a frictionless chain. 

Adding the small 0.5% correction for the uncompensated part of the friction on 
the table and the table edge, we find a total correction of 2.5%. The uncertainty 
of this correction is hard to estimate, but our gut feeling is that the correction is 
semi-realistic and that it is better to make it than not. We choose a relatively large 
uncertainty of 1.2% in order to allow for the other neglected effects mentioned in this 
appendix. This uncertainty is much greater than the uncertainty of 0.4% in the initial 
overhang position of the chain from the table edge. The latter uncertainty can thus 
be neglected. 

Appendix B. Chain tension, normal and frictional forces 

If the table is not smooth, two forces of kinetic friction act on the xyz-slider. They 
are (a) a horizontal force acting on the moving chain segment of length Sd = dj sin?/', 
and (b) a tangential force perpendicular to the normal force N at the table edge. 
These forces will be referred to in the following as the table friction and the edge 
friction, respectively. We shall show in this appendix that their effect on our measured 
acceleration is quite small. Nevertheless, their proper treatment requires a certain 
appreciation of the dynamics of our falling chains that is interesting in its own right. 

Consider first the simpler xy-slider whose vertically falling motion is illustrated in 
figure iBlT a). where the triangle shown can be of any size. The links on the horizontal 
side of the triangle shown in the figure as a dotted line are originally part of the 
horizontal and stationary chain segment at the table edge. These links then fall along 
the side of the triangle shown by the directed dashed line. Hence the chain appears to 
fall along the third side of the triangle shown as a solid line in the figure. The triangle 
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Figure Bl. Link and chain motions for (a) the icy-slider, and (b) the on-table 
chain segment of the xy^-slider. 



is thus an isosceles triangle with the apical angle 9 defined in section [3l and two equal 
base angles 

ai = \(-K-9). (B.l) 

The linear motions of links and chain occur everywhere on the falling chain if the 
apical angle 9 remains constant in time. 

Very much the same situation also appears on the table for the setback xyz-slidei. 
Here a short chain segment of length is sliding on the table with the links moving 
along the directed dashed side of the triangle shown in figure [BlT b). while the chain 
itself appears to move along the solid side of the triangle. The triangle has an apical 
angle ip and two equal base angles 

a 2 = i(7r-# (B.2) 

The table friction in the laboratory frame is thus opposite in direction to the directed 
dashed side of the triangle . 

We now turn to the friction at the table edge. Experience shows that the chain 
direction defined by tf> « 45° actually persists for some distance around and below 
the table edge before the overhang angle 9 as 70° finally takes over. Thus one may 
visualize figure IBlf b) to be folded around the table edge and to continue for some 
distance below it before the apical angle changes to 9. The plane containing the edge 
friction at the sharp table edge can be obtained by first placing figure [BlT b) on the 
table top and then rotating it about the table edge by an angle 4>n ~ 45° until its 
normal direction coincides with the direction of the normal force N. Edge friction is 
then opposite in direction to the directed dashed side of that rotated isosceles triangle 
whose lower base angle is located at the point where N acts. Just for the sake of 
notational transparency, we shall denote this new base angle a^. 

We are now in a position to decompose forces along different laboratory coordinate 
axes. For the z-components at the table edge of the xj/z-slider, we find 

N cos 4>n — HkN cos 4>n sin — pkpsdg sin ct2 

= T Z = pLZ = pals (B.3) 

for a chain whose center-of-mass position is (X, Y, Z) . This equation can be solved for 
the normal force: 

AT Pd9 (s . sma 2 \ . . 

N= T~n : \ -+A*fc— — T ( R4 

w OMpdg, (B.5) 

where the numerical value is calculated for pk = 0.13, ip = 45°, and a>3 = 
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As another example, the chain tension acting on the fixed chain end at the origin 
O of the xyz-slider is 

T = T y = pLY + HkN cos 03 + p, k pd s g cos ct2, (B-6) 

where the terms can be calculated in terms of s and the various parameters. The 
numerical result is T ~ Q.Hpsg. 
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